Anterior-posterior gene expression differences in three Lake Malawi cichlid fishes with variation in body stripe orientation

Morphological differentiation among closely related species provides opportunities to study mechanisms shaping natural phenotypic variation. Here, we address variation in the orientation of melanin-colored body stripes in three cichlid species of the tribe Haplochromini. Melanochromis auratus displays a common pattern of dark, straight horizontal body stripes, whereas in Aristochromis christyi and Buccochromis rhoadesii, oblique stripes extend from the anterior dorsal to the posterior mid-lateral trunk. We first validated a stably reference gene, and then, investigated the chromatophore distribution in the skin by assessing the expression levels of the iridophore and melanophore marker genes, ltk and slc24a5, respectively, as well as pmel, a melanophore pigmentation marker gene. We found anterior-posterior differences in the expression levels of the three genes in the oblique-striped species. The higher anterior expression of ltk, indicates increased iridophore density in the anterior region, i.e., uneven horizontal distribution of iridophores, which coincides with the anterior dorsalization of melanophore stripe in these species. The obliqueness of the horizontal body stripes might be a result of distinct migratory or patterning abilities of melanophores in anterior and posterior stripe regions which could be reflected by variation in the expression of genes involved in melanophore patterning. To address this, we investigated anterior-posterior expression levels of a primary set of candidate target genes with known functions in melanophore migration and stripe patterning in the adult zebrafish, and their related gene regulatory network. Among these genes, those with differences in anterior-posterior expression showed only species-specific differential expression, e.g., sdf1a, col14a1a, ifitm5, and agpat3, with the exception of fbxw4/hagoromo (differentially expressed in an oblique-and the straight-striped species). In summary, distinct anterior-posterior gradients in iridophore density found to be more similar characteristic between the two oblique-striped species. Furthermore, the species-specific differential expression of genes involved in stripe patterning might also implicate distinct molecular processes underlying the obliqueness of body stripe in two closely related cichlid species.